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Abstract 

We discuss the quantum paramagnetic phases of Heisenberg antiferromag- 
nets on the 1/5-depleted square lattice found in CaViOg. The possible 
phases of the quantum dimer model on this lattice are obtained by a map- 
ping to a quantum-mechanical height model. In addition to the "decoupled" 
phases found earlier, we find a possible intermediate spin-Peierls phase with 
spontaneously-broken lattice symmetry. Experimental signatures of the dif- 
ferent quantum paramagnetic phases are discussed. 
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The recent observation of a two-dimensional gapped quantum paramagnet in CaV40g 
has stimulated a number of theoretical studies IQ-Q of the spin 5=1/2 antiferromagnetic 
Heisenberg model with Hamiltonian H = '}2<ij> JijSi ■ Sj, and Si a spin S operator on 
each site i of the 1/5-diluted square lattice (Fig 1) found in this material. Of particular 
interest is the nearest-neighbor model, which has two exchange constants Ja, Jbc (see 
Fig 1). Theoretically, one possibility for the ground state is an ordinary Neel state in 
which the spins are ordered in opposite directions on the two sublattices. Alternatively, 
for Ja ^ Jbc or Jbc ^ Ja, the ground state is paramagnetic and its excitations are 
well understood; the lattice can be divided into essentially decoupled pairs {Jbc ^ Ja) 
or quadruplets {J a ^ Jbc) of spins. However, the real material is near neither limit [Q, 
and if the ground state is not Neel ordered, then there are a number of experimentally 
important questions about the resulting quantum paramagnet state. Among them are: does 
the ground state always have the full symmetry of the underlying lattice, as do the two 
states mentioned above, or can this symmetry be spontaneously broken? If the latter is the 
case, what is the nature of the excitations above the gap — does the system have unbound 
spin-1/2 spinous, or are they permanently bound into excitations with integer spin? Here, 
we will address these questions using the framework of earlier work on quantum paramagnet 
states of two-dimensional antiferromagnets ^-0]. We will find that while some of the physics 
is similar to that studied earlier on other bipartite lattices (the square and the honeycomb), 
the 1/5-diluted square lattice also displays some interesting new phenomena. 

Earlier work considered both unfrustrated and frustrated 0] Heisenberg antiferro- 
magnetic Hamiltonians on the square and honeycomb lattices. Based on various methods, 
a general theory of the quantum paramagnet phases of such systems was obtained, and also 
applied to the triangular and kagome lattices |p. The results are: (i) on bipartite lattices, 
with an unfrustrated Hamiltonian, the spinous are always confined to form integer-spin exci- 
tations, and the ground state must exhibit some form of spin-Peierls ordering that breaks the 
symmetry of the lattice, except possibly in the case when 25* is a multiple of the coordination 
number z of the lattice, when a nondegenerate state with the full lattice symmetry that has 
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only integer-spin excitations is possible 0; (ii) for frustrated systems, such as those with 
next-neighbor couplings on bipartite lattices or on non-bipartite lattices P, a ground 
state without spin-Peierls order can occur regardless of the value of 2S (mod z), with un- 
confined spinon excitations, which behave as bosons [0. All these states have an energy 
gap for local excitations, except at the second-order zero-temperature transitions that occur 
between some of them. 

The earlier results for bipartite lattices use the equivalence under lattice symmetries of all 
the links (and thus, nearest-neighbor interactions) of the lattices concerned. This property 
does not hold for the 1/5-diluted square lattice, and for the spin-1/2 antiferromagnet there 
are hmiting cases Ja ^ Jbc, and Jbc ^ Ja, in which it is evident that the ground state is 
invariant under all the symmetries of the lattice. However, these states fit into the general 
picture if we remark that, for non-Bravais lattices, it is clearly natural to consider the spin 
per unit cell as the analog of the spin per site on a Bravais lattice. For the present case, 
our lattice has a unit cell of 4 sites (unlike the statement by Ueda et al. 0) and we have 
possible total spin values of 0, 1, or 2 on a lattice with square symmetry, so a nondegenerate 
ground state along the hues of Ref. [^] is possible for the spin 2 case; if the sites of the square 
lattice are taken to be the A plaquettes, this can be interpreted as the ground state of the 
decoupled pair limit {Jbc ^ J a)- The decoupled quadruplet limit {J a ^ Jbc) is obviously 
the case of spin on each A plaquette. On the other hand, when the interactions are more 
nearly equal, a full analysis that treats each site as distinct may be more appropriate. We 
will show here that the limiting results, and a spin-Peierls phase, emerge in a unified way 
by a careful use of our earlier methods. 

Among the various approaches ||5|J^,pll-|T3| to the unfrustrated case that lead to the 
same effective action (all of which can be applied here), one of the most appealing is the 
quantum dimer (QD) model p!^^. This model represents the low-energy spin-singlet states 



by forming singlet "valence bonds" out of two spins of 1/2 on nearest-neighbor sites. Spin 
S* at a site can be represented by the symmetrized product of 25* spins of 1/2, and hence 
there must be 25* bonds ending at each site. Here, however, we will consider S* = 1/2 (the 



properties for other S are quite similar) and in this case each basis state of the QD model 
is associated with a close-packed dimer covering of the lattice under consideration, each 
dimer representing a valence bond. The effective Hamiltonian in this space consists of local 
diagonal potential energy terms, and off-diagonal terms which produce local rearrangements 
of the dimer packing. 

It has been shown [|l^,^|l5l that the imaginary time path integral of the QD model can 
be recast as a "height" or "roughening" model in space and imaginary time. The mapping 
relies on a one-to-one mapping between every state of the QD model and the equivalence 
classes of a configuration of heights, ha (on the sites, a, of the dual lattice (Fig 1)), under 
ha ^ ha + I. For the 1/5-depleted square lattice we have 

ha = na + (a (l) 

where is an integer which fluctuates from site to site, while (a is a fixed fractional offset: 
(a = 0,a,—a for a G A,B,C respectively, with 1/4 < a < 1/2 a fixed real constant. All 
configurations of ha which satisfy the constraint 

l^a ^ < 1 for every nearest-neighbor pair a, b (2) 

are permitted (two sites are nearest neighbors if their plaquettes on the direct lattice share 
a link). To map the ha configurations onto states of the QD model, examine the value of 
ha — hh for every pair of nearest neighbors, and if \ha — hi,\ > 1/2, only then place a dimer 
on the link of the direct lattice shared by the plaquettes around a and b. It can be shown 
that our choices for the offsets (a ensure that there is exactly one dimer terminating at each 
site of the direct lattice. 

It is helpful to examine the heights associated with some regular dimer coverings of 
the 1/5-depleted square lattice. The covering (I) ("fiat") in Fig 2 corresponds to /i^ = 0, 
hs = a, and he = — «; notice that this is the unique covering which is invariant under all 
the symmetries of the underlying lattice. The covering (II) ( "spin-Peierls" ) corresponds to 
hji = 0, hs = <y, and he = 1 — a; this covering has a partner, under the symmetry operations 
of the lattice that interchange B and C, which has hA = ^, hs = <y, he = 1 — a. 
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The path integral of the quantum mechanical height model involves summing over 
spacetime-dependent configurations of the heights, subject to the constraints, with an action 
related to the effective Hamiltonian. We will next write down a phenomenological effective 
action, S, that correctly describes the height model at long length and time scales; the same 
method can also be used for the classical dimer packing problem simply by making all the 
fields time-independent. We begin by using Poisson summation to promote the field ha, 
which can only take the discrete set of values (P, to a field Xa which can take all real values 
(the (a offsets in (|l]) require only a slight modification of the usual method [|6i|l2|). The 
constraints (|]) are "softened" by adding appropriate potential energy terms in the action; 
we will argue later that this softening cannot modify the main qualitative features of the 
results. We have: 



S= dr 



J2 ^abiXa - Xbf + J2 {^™ (drXa)'^ " 2/a COs(27r (Xa - Ca))} + • • • 



(3) 



_<ab> 

The terms shown are only representative; any term obeying the "translational" symmetry 
of the height model, Xa Xa + 1, is permitted. The coupling constants, Kat, K^a and 
Ha, depend only on the sublattice label of the sites, a, b, and the symmetry of the lattice 
requires that Kab = Kac, Vc = Vb etc. 

We now change variables from the three sublattice fields xa, Xb, Xc to the fields 
X2, X3 which are related by xa = Xi + X2, Xb = Xi - X2 + Xs, Xc = Xi - X2 - Xs- Now 
the translational symmetry of the height model affects only X15 Xi ^ Xi + 1; aiid this will 
lead to important simplifications. It is useful to note here the values of Xi, X2, Xs for the 
regular dimer coverings considered earlier. The state (I) ("fiat") in Fig 2 now corresponds 
to xi = 0, X2 = 0, Xs = ct- The two spin-Peierls states ((II)) have xi = V^j X2 = —1/4, 
X3 = « — 1/2, and Xi = 3/4, X2 = 1/4, Xs = ~ 1/2- These values suggest that X2 plays 
the role of the spin-Peierls order parameter, and that a non-zero mean- value of X2 implies a 
spontaneous breaking of the B ^ C interchange symmetry. 

The invariance of X2 and Xs under translations in height space implies that "mass" terms 
like ^2X2 + ^3X3 can, and do, appear in the action S. It is therefore safe to integrate out 



X2 and xs, which obey 

X2 ~ -sm(27rxi) , Xs ~ cos(27rxi), (4) 

at their saddle points, and obtain an effective action for the single scalar field xi- We take 
the spatial continuum limit for x\ and obtain 

Si = j drd^x [K{Wxif + Kr{drX2y - Vi cos(27rxi) - 2/2 cos(47rxi) - • • •] (5) 

Additional terms with more gradients of Xi, or cosines of higher integral multiples of 27rxi 
are also present. Note that, while, for the square (honeycomb) lattices considered earlier 0, 
lattice symmetries require that the effective potential for xi contain only cosines of integral 
multiples of Snxi (67rxi) (a fact intimately linked to the ubiquity of spin-Peierls order), 
here, lattice symmetries only impose the requirement that the effective potential be even in 
Xi, and we will see that phases without spin-Peierls order are also possible. 

As Si describes a three-dimensional interface, the interface must always be fiat, and the 
symmetry Xi — ^ Xi + 1 is spontaneously broken, even if all yi, 1/2, ■ ■ ■ , are zero. The value of 
(Xi), along with (Q), identify the state of the interface, and also of the QD model. A simple 
picture of the phases and phase transitions is obtained by truncating the on-site periodic 
potential to the two cosine terms explicitly displayed in (1), and minimizing the energy; 
more general potentials have the same qualitative features. The results are shown in Fig 2. 
There are three possible phases (modulo global translations of Xi by integers): 

(I) (xi) = 0: This is the "fiat" phase. There is no broken lattice symmetry, which is 
consistent with (X2) = 0. 

(II) < (xi) < 1/2: This state has spin-Peierls order and has spontaneously broken the 
B ^ C interchange symmetry because (X2) 7^ 0. The partner state is obtained by xi ~^ 

1 - Xl, X2 -X2, X3 X3- 

(III) (xi) = 1/2: this phase is similar to the "disordered flat" phase found in earlier work on 
two-dimensional interfaces |]TB|], so we will use that terminology [JT^. The interface is fiat on 



large scales, but each Ha fluctuates between two neighboring values with equal probability for 
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each, and correlations between /i^'s at different sites decay exponentially with separation. 
In the QD language, each A plaquette has two dimers which resonate between the two 
possible orientations, with correlations between the orientation of two plaquettes decaying 
exponentially. There is again no broken symmetry as (X2) = 0. Both the phases (I) and (III) 
are invariant under all lattice symmetries. They are nevertheless distinct states |0 which 
cannot be continuously connected. There is a non-trivial "topological order" , measured by 
the mean height, which distinguishes them. 

It is clear that phases (I) and (III) correspond to the Jbc ^ Ja and Ja ^ Jbc limits 
of the underlying antiferromagnet, respectively. Upon interpolating between these limits, 
we move along a section through the phase diagram, and there are two basic possibilities: 
(a) There is a direct first order transition between phases (I) and (III) as occurs in Fig 2 
for 1/2 > 0. (b) The spin Peierls phase ((II)) appears in between (I) and (III), as is the case 
in Fig 2 for ?/2 < 0. A third possibility is, (c), an intermediate phase with Neel long range 
order. This phase clearly lies beyond the scope of the QD model, and we expect that it can 
undergo direct second-order transitions to any of the three quantum paramagnetic phases, 
though if Neel and spin-Peierls phases are both present, we expect them to be adjacent. 
Finally, we also remark that the paramagnetic and Neel phases are stable under addition 
of not too large, non-nearest-neighbor, frustrating exchange interactions, although the Neel 
region is expected to shrink in size as frustration increases, as in Ref. 0. 

Beyond mean- field theory, fluctuations are expected to be relatively innocuous. First, 
as the interface is fiat in each phase, the softening of the constraint (|]) is not expected to 
have serious consequences; indeed, imposing (^ rigidly can only make the interface flatter. 
A direct transition between (I) and (III) (or other fiat phases with different mean heights) 
must always be first-order when the dimension of spacetime is greater than 2, because of the 
spontaneous breakdown of translational symmetry in xi even when the cosines are absent. 
In the spin-Peierls phase, the mean height varies continuously, and we expect a = 3 Ising 
transition to both (I) and (III), though first order behavior is not ruled out. The Ising-like 
behavior can best be understood by expanding the action in powers of xi about zero (or 
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1/2) to obtain a Landau- Ginzburg- Wilson 0^ action, with xi (or Xi ~ 1/2) playing the role 
of 0; the periodicity can be neglected here, since there is a much larger energy barrier for 
fluctuations changing x by ±1. 

To describe excitations of the antiferromagnet with non-zero spin, we use earlier re- 
sults Q, which go beyond the QD model, and show that a flat interface implies that spinous 
are confined. For the present model, this means that all three paramagnetic phases have 
only integer-spin excitations. 

We now turn to the classical, finite temperature (T) dimer model, whose height repre- 
sentation is the action Si without the time variable. The phase diagram in {K, yi, ?/2)-space 
is quite complex, but can be understood by fairly standard methods, so we will just sum- 
marise the results. The variable K ^ 1/T. At low T, a constant T section of the phase 
diagram near yi, y2 = resembles Fig. 2; the (I)-(III) transition is still first-order, while 
the (I)-(II) and (II)-(III) transitions are expected to be two-dimensional Ising transitions. 
For T greater than some critical value, the first-order transition line becomes second-order 
in the region adjoining the point where the three phase boundaries meet. There are then 
two Kosterlitz-Thouless (KT) critical points on the (I)-(III) boundary, one where the be- 
havior turns from first- to second-order, the other where the second-order portion meets 
the two Ising phase boundaries. The second-order portion of the (I)-(III) boundary is a 
region of continuously- varying exponents; the point 2/1 = 1/2 = lies in this interval . As 
T is raised further, a "rough" phase, with logarithmic correlations of Xi, appears around 
2/1 = 2/2 = 0, separated from (I) and (III) by a standard roughening transition, which is 
again of KT type [|19|. While completing this paper, a paper appeared by Weichman and 



Prasad |2^, who analyzed the two-dimensional height model with the two-cosine potential, 
finding similar results, but applied to roughening (see the inset to their Fig. 4). In the 
literature on classical two-dimensional roughening, a line of continuously-varying exponents 



is known as a preroughening transition |T^, and the phase II as 6'-disordered-flat p0|j21[] . 

For the antiferromagnetic Hamiltonian we started with, the situation is complicated by 
thermal excitation of nonzero spin excitations, and the finite T behavior remains unclear to 
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us, though we are confident that the three phases are separated by phase boundaries at low 
T, and that the lattice symmetry is broken below a finite in the spin-Peierls region. 

An experimental signature of the spin-Peierls phase (II) would be an accompanying lat- 
tice distortion, in which the links on which the spin-spin correlation is strongest, as shown 
in Fig. 2, would be shortened relative to the others, reducing the symmetry of the lattice. 
This would set in below the transition temperature Tc set by the antiferromagnetic inter- 
actions. It is not yet clear to us how the other phases (I) and (III) can be experimentally 
distinguished from one another, although in principle a jump in the lattice constants should 
be observable when the first-order transition at low T is crossed. In all three phases, the 
elementary spin-carrying excitations have integer spin, and spin 1 is most likely. This should 
have a clear signature in polarized neutron scattering: the dynamic structure factor S{k, uj) 
{k and lu are the transferred momentum and frequency) should have a quasiparticle delta 
function ~ 5(a; — e{k)) at T — 0, where e{k) is the dispersion relation of the spin 1 quasipar- 
ticle. This signature is analogous to those of 5" = 1 Haldane gap antiferromagnets in d — 1, 
and is common to all phases in which spinous are confined, of which a survey was given at 
the beginning of this paper. 

To conclude, our main prediction is the appearance of the spin-Peierls phase, as shown 
in the phase diagram in Fig. 2. 

We thank R.R.P. Singh and C.L. Henley for stimulating our interest on this problem, 
and especially the latter for pointing out that height formulations of the QD model could 
be extended to the 1/5-dcpleted square lattice. This research was supported by NSF grants 
nos. DMR-92-24290 and DMR-91-57484. 
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FIG. 1. The 1/5-depleted square lattice, and the three sublattices {A, B, C) of its dual lattice. 
The exchange J a acts between spins on a link in a A plaquette, while Jbc acts on links shared by 
the B and C plaquettes. For Jbc ^ Ja, the ground state is in the "flat" phase of Fig. 2, and for 
Ja S> Jbc it is in the "disordered flat" phase. 
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FIG. 2. Mean-field phase diagram of the model defined by 5*1 (Eqn (^)) with the potential 
truncated to two cosines. The thick line is a first order transition, while the thin lines are second 
order. Effects of fluctuations are discussed in the text. 
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